B/C ratio and the PAMELA positron excess 
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We calculate the diffuse intensity of cosmic ray (CR) nuclei and their secondaries in the Boron- 
Carbon group produced by supernova remnants (SNR). The trajectories of charged particles in the 
SNR are modeled as a random walk in the test particle approximation. Secondary production by 
CRs colliding with gas in the SNR is included as a Monte Carlo process, while we use Galprop to 
account for the propagation and interactions of CRs in the Galaxy. In the vicinity of a source, we 
find an approximately constant B/C ratio as a function of energy. As a result, the B/C ratio at 
Earth does not rise with energy, but flattens instead in the high energy limit. This prediction can 
be soon tested by the AMS-2 collaboration. 
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I. INTRODUCTION 

Measurements of the nuclear composition and of the 
antimatter fraction of cosmic rays (CR) are one of the 
main tools to understand their origin Because the 
propagation of charged CRs in the turbulent component 
of the Galactic magnetic field is diffusive up to energies 
E/Z < 10 17 eV, properties as the average injection spec- 
trum of CRs or the diffusion tensor and its energy de- 
pendence have to be determined indirectly. The qual- 
ity of these measurements has been improved in the last 
years considerably, for instance by the determinations of 
the intensity of CR protons and He performed by the 
PAMELA collaboration @]. In particular, the PAMELA 
collaboration found a slope change of the CR spectrum, 
which is difficult to explain in the simplest models Q. 
Similar deviations from a pure power-law even below the 
knee region were found using gamma-ray observations of 
selected molecular clouds in the Gould belt d, Q . 

Previously, the PAMELA collaboration measured the 
positron fraction in CRs and found a rapid rise from 10 
to 100 GeV @. Fermi LAT measurements confirmed this 
behavior in the energy range 20-200 GeV Q- In con- 
trast, the antiproton ratio measured by PAMELA de- 
clines above 10 GeV [8|, consistent with expectations. 
Such a rise of the positron fraction is most naturally ex- 
plained by the injection of high-energy positrons by local 
sources as nearby pulsars Q: Since electrons loose fast 



energy, the high-energy part of the 



e + spectrum 



is dominated by local sources. Moreover, electromag- 
netic pair cascades in pulsars result naturally in a large 
positron fraction together with a "standard" antiproton 
flux. 

In this work we study an alternative explanation for 
the rising secondary fraction which was put forward orig- 
inally for positrons in Ref. [lOj: Secondaries created by 
hadronic interactions in the shock vicinity participate in 
the acceleration process and therefore it was suggested 
that they have a flatter energy spectrum than primary 
CRs. In Ref. [lOj, it was estimated that the resulting 
positron fraction can explain the PAMELA excess and 



increases up to 50%, while subsequently this mechanism 
was applied to antiprotons in Ref. ll| and to the B/C 



and Ti/Fe ratios in Refs. [12 1 



The present work examines the acceleration and pro- 
duction of nuclei in SNRs. Following closely in spirit our 
previous studies [H, [13], we use a Monte Carlo (MC) 
approach calculating the trajectory of each particle indi- 
vidually in a random walk picture. This makes it easy to 
include interactions and the production of secondaries. 
As our main result, we predict a flat B/C ratio in the 
source vicinity, without a rise at high energies. At Earth, 
the additional contribution from SNRs to the secondary 
fluxes only flattens the B/C ratio but does not lead to a 
rise in the high energy limit. Consequently, an increase 
in the antiproton-proton ratio cannot be explained by as- 
trophysical models and could be used as a signature of 
DM. This prediction can be tested soon by the AMS-2 
collaboration 11511 . 



II. SIMULATION PROCEDURE 

As we are follo wing closely the approach presented al- 
ready in Refs. [l3|, Il4| . we recall only briefly our method- 
ology. We model trajectories by a random walk in three 
dimensions with step size Iq(E) determined by an energy- 
dependent, isotropic diffusion coefficient D. We assume 
as usual in the test particle approach that the diffusion 
of a nucleus with charge Ze and momentum p proceeds 
in the Bohm regime, 
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where B denotes the turbulent magnetic field. As a draw- 
back of our Monte Carlo approach, the step-size in the 
random walk decreases as Iq cx 1/Z, requiring larger com- 
puting time for nuclei with higher charges Ze. Therefore 
we restrict our analysis to relatively light nuclei, concen- 
trating on the B/C ratio. 

We do not consider any feedback of CRs on the shock 
or the magnetic field. As the strength of the magnetic 
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field is assumed to change from a phase of amplification 
to damping, we use two models for the time evolution of 
the magnetic field B: In the first case, we use a constant 
turbulent magnetic field with B — 1 fj,G during all the 
evolution of the SNR. As the second choice, we assume 
strong field amplification in the early phase, with B = 
100 fiG before the transition to the Sedov- Taylor phase 
at i* = 240yr, and subsequent damping to B = 1/20 [iG 
at t ^ . 

For the position r s h and the velocity u s h of the SNR 
shock we use the n = case of the analytical solutions 
derived in Ref. [l6(. These solutions connect smoothly 
the ejecta-dominated phase with free expansion r s h oc t 
and the Sedov- Taylor stage r s h oc t 2 / 5 . The acceleration 
of CRs is assumed to cease after the transition to the 
radiative phase at the time i max . We use an age-limited 
scenario for the CR escape where most CRs are accumu- 
lated downstream and are released at i max . The SNR 
is modeled by the following parameters: We choose the 
injected mass as M j = 4M , the mechanical explosion 
energy as _E snr = 5 x 10 51 erg, and the density of the ISM 
as riisM = 2 cm -3 . The end of the Sedov- Taylor phase 
follows then as £ max = 13.000 yr. We model the injection 
rate proportional to the CR pressure [TtJ , 



N ex R 2 v° h 6(E - E Q )S{r - r sh ) 



(2) 



with a = 3. We choose the isotopes according to the 
source abundances of Ref. |l|. Lifetimes for radioactive 
isotopes are taken from the database [Hj]. We use for 
the simulation of nuclear reactions the cross sections and 
decay tables extracted from GALPROP which in turn 
are based on Refs. fl9l I20I] . 

We describe the propagation of cosmic r ays in the 
Galaxy with the help of the Galprop code 0, [2l[ . More 
specifically, we consider plain diffusion, neglecting re- 
acceleration process and fixing all relevant parameters, 
like the normalization and the energy-slope of the diffu- 
sion coefficient, from a fit to HEAO data [22| to both 
primary and secondary nuclei in the 1 -j- 30 GeV range. 



III. NUMERICAL RESULTS 

In Fig. [TJ we show the energy dependence of the B 
( 10 B+ n B) and C ( 12 C+ 13 C) spectra E 2 dN/dE at i max 
for the two cases of magnetic field evolution. For a con- 
stant magnetic field, the spectra dN / dE of both primary 
and secondary CRs are well approximated by an E~ 2 
shape 1 . The corresponding B/C ratio shown in Fig. [5] by 
the solid line remains constant up to the acceleration cut- 
off: Because the energy per nucleon E/A in a spallation 



1 The small "bump" near the acceleration cutoff -E m ax ~ few X 
10 13 ZcV is formed by the small overdensity of the most ener geti c 
CRs in the acceleration zone, as discussed in more detail : 
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FIG. 1: Boron (dashed blue) and Carbon (solid red) energy 
spectra E 2 dN/dE at the source as a function of energy per 
nucleon for two different magnetic field models: left panel 
constant and right panel time-dependent magnetic field. 
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FIG. 2: The B/C ratio at the source as a function of energy 
per nucleon for the constant (solid red) and time-dependent 
(dashed blue) magnetic field case. 



reaction is approximately conserved and the simple rela- 
tion E/Z ~ 2E/A holds, the rigidity E/Z of the reaction 
products stays approximately unchanged. Hence, the ac- 
celeration process is not influenced by the conversion of, 
say, a 12 C into a 10 B nuclei: The secondary nucleus fol- 
lows the same trajectory in a given magnetic field con- 
figuration as its parent would do, if no conversion took 
place. As a consequence, the boron and carbon spectra 
as function of energy per nucleon have practically identi- 
cal shapes, with their relative normalization determined 
approximately by the interaction depth in the SNR: The 
latter is given by r = ct max Rn lsm ac ~ 0.01 assuming 
ere ~ 50mbarn as typical value for the fragmentation 
cross section of a carbon nuclei. 

In the second model for the magnetic field evolution, 
the spectral shape is more complicated: Nuclei that were 
injected early are accelerated up to few xl0 15 ZeV, while 
the bulk of CRs injected when the turbulent magnetic 
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FIG. 3: Diffuse intensity of Boron (red dashed) and Carbon (red solid) at the position of the Earth for the constant (left 
panel) and time-dependent (right panel) magnetic field case. Additionally, in the left panel the B and C intensity including 
only secondary production in the ISM are shown by black dot-dashed and black doted lines respectively. Experimental data: 
HEAO 22] - filled circles, CREAM [H - open squares, TRACER - filled squares. 
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FIG. 4: The B/C ratio as a function of energy per nucleon 
at the position of the Earth for the constant (red solid) and 
time-dependent (blue dashed) magnetic field case; the green 
dot-dashed line shows the B/C ratio including only secondary 
production in the ISM . Experimental data: HEAO [22| - filled 
circles, ATIC 12511 - open stars, CREAM [2g] - open squares, 
and TRACER [2J] - filled squares. 



field is damped has a cutoff around 10 12 ZeV. However, 
due to the conservation of the rigidity by the spallation 
process, the energy spectra of B and C are remarkably 
similar: The B/C ratio is again energy-independent and 
coincides with the one obtained for the case of a con- 
stant magnetic field — as it should, since the interaction 
depth is the same in the two cases. Thus, our conclusion 
on the energy-independence of the primary to secondary 
CR ratios is of general character, being independent of a 



particular scenario for the source magnetic field. 

As the next step, we calculate the CR propagation 
in the Galaxy using the Galprop code, and compare in 
Fig. [3] the resulting B and C intensity at Earth with ex- 
perimental data. For the spatial diffusion coefficient in 
the Galaxy, we used D XX {E) = D (E/10 GeV) 065 , with 
D = 4.8 x 10 28 cm 2 /s. In the left panel of Fig. H we 
show additionally the B and C intensity including only 
secondary production in the interstellar medium (ISM), 
i.e. neglecting the SNR contribution to the B intensity. 
The SNR contribution to the Boron intensity remains 
small in the whole energy range considered, independent 
of the models used for the magnetic field evolution in 
the SNR. At Earth, the obtained B/C ratio shown in 
Fig. |4] coincide with the pure propagation results up to 
~ few x 100 GeV/nucleon; at higher energies the ratios 
slowly flatten to the per cent level, i.e. to the respective 
source value. 



IV. CONCLUSIONS 

We calculated the energy spectra of nuclei in the B-C 
group and their secondaries produced in a supernova rem- 
nant (SNR) using a simple random walk picture. In con- 
trast to a previous prediction [l2[ that the B/C ratio rises 
sharply above energies ~ TeV/nucleon, we found that the 
ratio remains approximately constant in the vicinity of 
a source. Our results suggest that the contribution of 
SNRs to the intensity of secondary CRs is subdominant 
up to very high energies. 

In our simulations we used an age-limited escape of 
CRs. Implementing a free-escape boundary may change 
our results, except if the boundary is placed at the dis- 
tance oc E. In particular, assuming the same position of 
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the free-escape boundary for all e nerg ies will lead to a 
softening of the secondary spectra [27| . lowering further 
the B/C ratio. 

This prediction can be soon tested by the AMS-2 col- 
laboration The absence of a rise in the B/C and 
Ti/Fe ratio would support the arguments of Refs. [l3l.[l4T| 
that reacceleration does not lead to a change in the sec- 
ondary to primary ratios. As result, a rise of the an- 
tiproton ratio could be used as signature for dark matter 
searches. 
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